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ABSTRACT 


The distribution of gas satisfying hydrostatic and Poisson conditions 
with distance above the galactic plane is derived and compared with 
Schmidt's observations of Pg(z). Equipartition magnetic (B z /8rr) and 
cosmic-ray (Pc-r) components are assumed. Pg( z ) is calculated for two 

limiting cases: (i) the total mass density at the galactic plane 
p go + p* 0 is equal to the sum of the observed gas and star densities 
in the solar neighborhood (0.089 t^/pc^), and (ii) Oort's data on the z 

component of the galactic gravitational acceleration K z are used, in 

3 

which case we implicitly assume that Pg 0 + P*o := 0.15 Mg/pc . Using 


recently observed values for B 0 , P c _ ro , and p g0 , the rms z turbulent 


2 1/2 

gas velocity dispersion < v tz > is found. 


An expression is derived for the half- thickness ^\j r i of the 

1/2 2 2 

equilibrium disk of galactic gas: z^/2 * Q/(Pgo + P*o) where Q - < v tz > + 
2 

B 0 /8rrpg 0 + Pc-ro/Pgo* Using the Innanen galactic mass model, the 

half-thickness is computed as a function of distance from the galactic 

center over the range 4kpc^RS14kpc and compared to the observations of 

McGee and Milton. The observed increase in half- thickness at distances 

beyond the solar distance R Q is reproduced theoretically. A galactic 

mass model is derived, using the observed layer half-width and a constant 

value for Q. Comparison of the two mass models indicates that, because 

z\/2 is observed to remain essentially constant in the region 

, 1/2 

4kpc£R*!R 0 , Q increases like (p gQ +p* 0 ) 


with decreasing R. 
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I„ INTRODUCTION 

In recent years observations of interstellar molecular absorption 
and emission lines (e.g. OH, H 2 O, NH 3 , CH, CN, CO etc.), interstellar 

extinction, and 21 -cm radiation have aided in determining the distribution 
of the gaseous component in our galaxy. This paper discusses, among 
other topics, two observational findings relating to the equilibrium 
state of the galactic gas disk: (i) the run of gas density Pg with 
distance above and below the galactic plane (Schmidt 1956), and (ii) 
the rapid increase in half-thickness of the gas layer at distances 
from the galactic center greater than the solar distance (McGee and 
Milton 1964), and its relatively uniform thickness at distances 
between about 4 kpc and 10 kpc (Schmidt 1956; McGee and Milton 1964; 

Kerr 1969). To compute the equilibrium state of the gas disk, the 
Poisson equation and the hydrostatic equilibrium equations for the 
gaseous and stellar components of a two-fluid, static, plane-parallel 
mixture are used. Each component is assumed to be isothermal but 
each has its own characteristic 'temperature'. The analysis includes 
one-dimensional magnetic and cosmic- ray components, each constrained 
to vary with p g (z) according to equipartition of energy. 

II. THE EQUILIBRIUM DISTRIBUTION OF 
GAS ABOVE THE GALACTIC PLANE 

The equilbirum configuration of the galactic gaseous, stellar, 
magnetic, and cosmic-ray components may be pictured as follows. A 
plane-stratified (in z) infinitely-extended (in r) non-rotating gas 
slab is immersed in a star slab of similar description, the gas and 
star densities attaining maximum values at the mid-plane z=0. The 
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weight of the gas layer confines a one-dimensional magnetic field 
and cosmic-ray gas, the direction of the field assumed parallel to the 
plane z=0. Early optical polarization measurements (Hiltner 1956) gave 
support to such a simple magnetic topology. However, more recent optical 
polarization measurements (Mathewson 1968; Mathewson and Nicholls 1968) 
on large numbers of stars indicate the presence of a helical component 
with a pitch angle of about 7°. 

The Poisson and hydrostatic equilibrium equations expressed in 
cylindrical geometry may be written as follows: 


= -4rrG(p„+p*) (1) 

dz 2 S 

i(P g + B 2 /Su + P c - r ) - P g f| (2) 



= D dcp 
p *d^ 


(3) 


where p g and p* are the gas and star densities, P g , P*, P c _ r , and 

O 

B / 8nr are, respectively, the gas, star, cosmic-ray, and magnetic 
pressures, and cp is the total gravitational potential, contributed to 
by both gas and stars. The terms containing derivatives with respect 
to r in the expressions for V and V are omitted. Their contribution 
is estimated to be less than 10% of the total at the sun's distance 
from the galactic center. The equation of state relating the 
pressure and density of the isothermal gaseous component can be 
written as 


P g = <v tz > Pg 


where <v tz > is the 


( 4 ) 


mean square z turbulent and thermal gas velocity 
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dispersion. Following Parker (1966, 1970) we assume that <v tz > and 
the ratios of magnetic and cosmic-ray pressure to gas pressure are 
independent of z, allowing us to write that 
B 2 (z)/p g (z) = B o 2 /p g0 


^c-r p g( z ) ~ ^c-ro^Pgo* 

where the subscript o designates that the quantities are evaluated at 
the galactic plane (z=0). With the aid of equations (5), equation (4), 
and a similar equation for the stellar component 

P* = <v*2> P*, (6) 

the two hydrostatic equilbrium equations assume the form 

= to < 7 > 

2 d d<P 

<v *z > dz Anp * “ dz’ ( 8 ) 

2 2 

where Q 2 denotes the quantity <v tz > + B c /8rrp go + P c _ ro /p go . We have 

explicitly assumed here that both the gaseous and stellar components 
of the two-fluid mixture are isothermal, i.e. that <v 2 2 > and <v* z > 
are independent of z. 

It is not necessary to use the Poisson equation and the stellar 
hydrostatic equilibrium equation explicitly to determine the distribu- 
tion p g (z)/p g0 if the gravitational acceleration ckp/dz = K z is known 
as a function of z. Oort (1960) determined this acceleration by 
analyzing star counts and radial velocities of K giants. Oort found 
the density of all matter in the solar neighborhood to be 0. 15 i^/pc 3 , 

significantly greater than the combined observed density of stars 
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3 

and gas which is 0.089 tfe/pc (Luyten 1968; Weaver 1970). More recently 

c 

Woolley and Stewart (1967), noting that the velocity dispersion of A 
type stars increases near the galactic pole with increasing magnitude, 
concluded that the density of all matter in the solar vicinity is 

3 

O.llF^/pc , considerably closer to the observed mass density than found 
by Oort. 

The question of the true density in the solar neighborhood and 
hence of the component of the gravitational acceleration perpendicular 
to the galactic plane is not yet fully settled. We choose therefore 
to follow two separate courses in calculating Pg(z), keeping in mind 

that the 'correct' approach may lie between the following two limiting 
cases. (i) The Poisson equation and the two hydrostatic equilibrium 
equations are solved simultaneously for p g (z), assuming that the total 
mass density at the galactic plane is represented by the sum of the 

O 

observed gas and star densities in the solar neighborhood (0D89 Mg/pc ) . 
Since the bulk of the stellar mass in the solar vicinity resides in 
the main sequence G, K, and M stars (Oort 1960) and since each type 
has nearly the same velocity dispersion of 18 km/sec in the z direction 
(Wehlau 1957; Woolley 1958), a unique value of 18 km/sec for <v* z >^^ 
will be assumed for all stars. (ii) The gas hydrostatic equilibrium 
equation is solved for Pg(z) using Oort's K z , and thus implicitly 

3 

assuming that p go +P*o = 0.15^/pc . This procedure differs from 
Parker's (1966, 1970) work in three ways. First, the simplifying 
approximation K z (z) = <K Z > is not employed. Second, we take B 0 , 

2 i/2 

P c _ro» an< * Pgo as known from the observations and derive <v tz > ' 
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by equating the observed half-thickness of the gas disk with the 

2 1/2 

calculated half-thickness . Parker assumed that B 0 , P c _ ro , and < v t2 > 
are known and derived p gQ . It is the opinion of the present writer 

2 l/o 

that <v tz > ' * is the least well-known of the four quantities for the 
following reason. 21-cm profiles have so far given information on 
only the component of the turbulent gas velocity dispersion in the 
galactic plane (6-7 km/sec) since most of the neutral hydrogren is 
confined to a thin disk, with the sun near the mid-plane. There is 
essentially no data for <v tz > ' , although there may be some theoretical 
justification in taking = <V£J>^^. We do not 

choose to make such an assumption. Third, the distribution of gas 

above the galactic plane p g (z) is derived once the appropriate value 

- . 2 1/2 . , 

of <v tz > is found. 

We proceed with (i) by equating the left-hand sides of equations 

(7) and (8) and integrating, with the result that 

2 2 

P* = P*o(Pg/Pgo)^ * z (9) 


Combining this with the Poisson equation then gives 

9 2 

d 2 cp Q /<V*2> 

= -%G[p g +p* 0 (p g /p g0 ) ]. 


The gas hydrostatic equilibrium equation is differentiated and combined 
with equation (10): 

2 2/2 

4npg/p go = - ~ [p g +p* 0 (p g /p go ) Q VCzr> ]. (11) 

dz^ 

Equation (11) may be reduced to a dimensionless form by defining 
four dimensionless parameters 
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a = Pg(z)/p go 

8 = z/H g where H g = Q/ 0rtGp go ) 


y p*o^Pgo 


6 = Q 2 /<v| z >, 

in which case equation ( 11 ) becomes 
d 2 


d02 


iwt = - 1 / 2 (a + 7 a 6 ) , 


( 12 ) 


(13) 


Equation (13) is solved numerically for a ( 8 ) with three values of 6 
and hence of Q. Q is chosen to be 5.0, 7.7, and 10.0 km/sec: < v £ > 1 ^ 2 
is chosen to be 18 km/ sec as mentioned above. The observed gas and 
star densities p go = 1 H atom/cm = 0.025 %/pc (Weaver 1970) and 

P*o = 0.064 f^/pc^ (Iiiyten 1968) are used. The resultant curves of 
Pg( 2 )/Pg 0 are presented in Figure 1 together with a smooth curve 

through Schmidt's (1956) data, the latter corrected by the change in 
galactic distance scale R 0 (new)/R Q (old) = 10.0 kpc/8.2 kpc = 1.22. 


When z < 200 pc the solution with Q = 7.7 km/sec agrees closely 
with Schmidt's observations. When z > 200 pc the observed curve lies 
above the theoretical solution, the discrepancy increasing with 
increasing z. In fact, the curve with Q = 10.0 km/sec fits the 
observed curve at z = 400 pc, an increase of 2.3 km/sec. Some possible 
causes of this discrepancy will be discussed shortly. 

Since Q is composed of turbulent, thermal, magnetic and cosmic - 
ray components, i.e. 

Q “ <v tz> / Srtpgo "*■ Pc-ro/Pgo* 


(14) 
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2 1/2 

< v tz > can be determined by employing the Q found from Figure 1 and 
the observed values of B 0 2 /8 tt, P c _ ro> and p gQ . Using a cosmic-ray 

pressure of 0.50 x 10 dynes/cm (Parker 1970), a magnetic pressure of 
0.36x10 dynes/cm corresponding to a field strength of 3.0x10 ^gauss 
(Verschuur 1970; Manchester 1971), and p gQ = 0. 025 1^,/pc 3 , we find that 
<v tz >1/2 = 3 * 00 Km/sec. 


Recent evidence from RR Lyrae star counts near the galactic 
center (Plaut 1970) indicates that R 0 may be significantly less than 

10 kpc. If we choose R 0 = 8.2 kpc and repeat the above procedure, 

2 1/2 

<v tz > ' turns out to be 0 km/ sec, i.e. there is no room for z gas 
turbulence. This result can clearly be ruled out. 

The second approach consists in solving the gas hydrostatic 
equilibrium equation using Oort's K z . We recall equation (7): 


Q2 ^ np g /p g° - ^ ~ K z‘ 


(15) 


Integrating between the limits z=0 and z=z we can solve for p„ ( z )/ p „ 0 ; 

9P Z 8 8 

= e 1/Q i 


Pg (2)/ Pgo = 


(16) 


Note that neither p* Q nor <v* z > appear here. Rather, they are 
contained implicitly in K z . The computed distribution Pg (z)/Pg D is 
shown in Figure 2 together with a smooth curve through Schmidt's data, 
once again corrected for the change in galactic scale length. 

Q = 9.84 kn/sec gives the closest fit to the observed curve for 
z < 200 pc. This solution follows closely the simultaneous solution 
of the Poisson equation and the two hydrostatic equilibrium equations 
computed above (procedure (i)) with Q = 7.7 km/sec, as a comparison 
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of Figures 1 and 2 shows. The higher Q arising from the use of Oort's 
K z is directly related to the higher density of all matter found by 
Oort (0.15 vs. 0.089 M^/pc^) at the galactic plane. Using Q = 9.84 km/sec 

O 

with the above quoted values for /8n, P c _ ro , and Pgo> we find that 

<v tz> 1/2 = ^-82 km/sec. We should emphasize that < v^ z >^^ is a function 
of the value chosen for R 0 , since a change in R 0 is reflected in the 

observed half- thickness of the gas layer. A decrease in Rq causes a 

corresponding decrease in half-thickness , and a smaller gas turbulent 

velocity is necessary to support the gas layer against gravity. 

2 1/2 

Specifically, if R Q = 8.2 kpc rather than 10.0 kpc , <v tz > becomes 
3. 51 km/ sec. 

Let us comment briefly here as to the outcome of approaches (i) 

and (ii). Only (ii) using Oort's K z (implicitly assuming that p go + p*o 

0.15 %/pc^) and R 0 = 10 kpc yields £ n t ^ e neighborhood of 

7 km/sec. All other assumptions result in z turbulent gas velocity 

2 1/2 

dispersions at least a factor of two smaller. Indeed, if -Cv^^ 7 = 

2 1/2 2 1/2 2 1/2 
<v tx > = <V t y> ' as some authors have suggested, and if <v tx > and 

2 1/2 

<v ty > X equal 6-7 km/ sec as observed by Westerhout (1956) and others, 
we may conclude that the evidence presented above favors the Oort mass 
limit as being substantially correct and R 0 close to 10 kpc. Table 1 

summarizes the results of the present discussion, giving < ' v tz' > ^'^ 
as a function of p go +p ^ and R Q . 

Recent high-resolution 21-cm observations (Kerr 1969) indicate 
that a smaller value for the layer half-thickness may be appropriate 
(200 pc), which is in better agreement with the layer width found for 
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Cepheids and HII regions (Kerr 1964). It is clear from our discussion 
above that less than 3 km/sec would be needed to support the 

gas layer to the observed width, even if R o =10 kpc. 

The theoretical distribution Pg( z )/Pg 0 f a lls off too rapidly with 
distances above the galactic plane greater than about 200 pc using 
either approach. We consider five possible explanations. (a) The gas 
is not in hydrostatic equilibrium above about 200 pc from the mid-plane, 
(b) The equipartition assumptions that B^/pg^v^^ and P c _ r /p <v t z > are 
independent of z break down when z > 200 pc. (c) The observations 
at z > 200 pc are incorrect. (d) Schmidt's data for z > 200 pc probably 
refers to tangential points significantly closer to the galactic center 
than R c , where increased 21-cm line profile widths (Schmidt 1956) may y 
indicate an increase in <v t From this fact alone we would expect 
more gas at large z. However, this is not the whole story (see 
Section III), since Pg 0 +P, V0 and thus K z is also expected to increase 

with decreasing R, which would tend to reduce the amount of gas at 
large distances from the mid-plane. We will show below that these 
two efforts cancel, so that this explanation is probably not the 
correct one. (e) According to Field, Goldsmith and Habing (1969) a 
critical pressure P c exists such that if Pg(z)> P c , the neutral gas is 
essentially in a cool dense phase and if P (z) < P c the gas is in a hot 

O 

rarefied phase. In addition they showed that Pg Q is such that the gas 

near the mid-plane should predominantly be in the cool phase. There- 
fore, since dPg/dz < 0 from hydrostatic considerations, there will be 



12 


a certain distance z c above the mid-plane at which Pg = P c . When 

j zj < z c the gas will be predominantly in the cool dense phase; when 
|z| > z c the gas will be mainly in the hot phase. From our calculations 
Pg = 0.22 atoms/cm^ at z = 200 pc (assuming pg Q = 1.0 atoms/cm^), which 

3 

agrees closely with the p c found by Field et al. (1969) = 0.21 atoms/cm 

3 

by Spitzer and Scott (1969) = 0.21 atoms/cm , and by Hjellming, Gordon, 

and Gordon (1969) = 0.19 atoms/cnr*. The divergence between the observed 

and theoretically derived distributions Pg(z)/pg Q above z=200 pc may 

therefore be caused by the presence of the high temperature (8000°K) 

2 1/2 

phase with large <V tz > ' that should predominate above z=200 pc. It 

is of some interest to note here that Kepner (1970) observed gas in 

spiral arms extending 1-2 kpc from the galactic plane. We suggest 

that pg (z)/p go be derived from the more recent 21-cm data so that these 

points may be clarified. 

III. THE HALF-THICKNESS OF THE 
EQUILIBRIUM GAS DISK 

Schmidt (1956) observed the thickness of the neutral hydrogren 
layer at several galactic tangential points and found that between 
about 3 kpc from the galactic center and R 0 , the layer half-thickness 
is surprisingly constant, with a value of 220 pc (268 pc on the new 
distance scale with R o =10 kpc). More recently Kerr (1969) has found 
a somewhat smaller constant half-thickness of about 200 pc over the 
region 4 kpc S R S Rq. Westerhout's (1956) isodensity contours of 
neutral hydrogen indicate that the half-thickness of the gas layer 
increases substantially with increasing distance from the galactic 
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center beyond the sun's position Rq. McGee and Milton (1964) confirmed 
the constant half- thickness at distances between about 3 kpc and R Q 
(on the old distance scale with R c = 8.2 kpc). Beyond R 0 they found 
the half-thickness to increase rapidly, reaching more than 1000 pc 14 kpc 
from the galactic center. We shall attempt to explain these observations 
from a theoretical point of view. 

We first recall equation (13): 

d 2 R 

— Urt* = - 1/ 2 (o'+yar ) . (13) 

d$ 2 


6 is typically on the order of 0.2-0. 3 from our discussion above and 

can be set equal to zero with little resulting loss of accuracy. 

Further, a is equal to pg(z)/p g0 , and since the latter may be approxi- 
_8 2 

mated by e p with good accuracy, equation (13) becomes 


2 2 
- — ira » - l/2(e"^ +y). 
dp2 

-ft 2 


(17) 


_Q 4 O 

Approximating e p by 1-8^ and twice integrating, equation (17) becomes 
Ira « - l/4(Hy)8 2 . (18) 

Setting a = Pg< z )/p g0 = 1/2, equation (18) can be solved for P 2 * 
q / 2.773N1/2 /1QA 

The half-thickness , i.e. the distance between the points above and 
below the galactic plane where the gas density is one-half its value 
at the plane, expressed in dimensional form, is written Zjy2 = 2 Pl/2^g : 

1/2 n 


: l/2 


/ 1 o 386\ 

\ nG / 


(P + P ) 
go y *o 


172 * 


( 20 ) 
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Equation (20) reveals that the half-thickness is directly proportional 

2 2 2 

to Q, where we recall that Q = <v tz > + B 0 /8rTp go + P c -ro/pgo» and is 
inversely proportional to the square root of the total mass density at 
the galactic plane. 

Taking z^/2 as known from the observations of McGee and Milton 
(1964), we may use equation (20) (i) to compare the observed and 
predicted distributions zjy 2 (R), using a galactic mass model for p^ Q (R)= 

p go (R) + p* 0 (R) and choosing that value of Q which gives the best fit 

to the observations at R = R 0 (where Q is assumed to be independent 
of distance from the galactic center), (ii) to derive a galactic mass 
model assuming an appropriate value for Q independent of R, and (iii) 
to derive Q(R) on the basis of a galactic mass model. Schmidt (1956) 
has explored (iii) for R < R 0 using his galactic mass model. We shall 
pursue (i) and (ii) and comment briefly on (iii). 

To pursue (i) the quantity p^o = Pg 0 +P*o (the total density of 
matter at the galactic plane) must be known as a function of distance 
from the galactic center. Interstellar obscuration by dust severely 
limits our knowledge of p^ 0 (R) at distances beyond about 1 kpc from 
the sun. Theoretically determined densities must therefore be used. 

We have selected the Innanen (1966) galactic mass model with Rq = 10.0 kpc 
and Vq = 2.52 km/ sec (Table 2). Q is chosen to be 7.25, 10.25, and 
13.25 km/sec. The three curves for that result are compared 

with the observations in Figure 3. 

At distances less than about 8 kpc from the galactic center, all 
three of the theoretically derived curves fall below the mean curve 
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through the observed points, the discrepancy increasing with decreasing 
R. This suggests that either a larger value is required for Q or that 
the total mass density at the galactic plane p^ 0 has been overestimated. 
The curve with Q = 10.25 km/sec exhibits the best agreement with the 
observed half- thickness of 268 pc at the sun's position (assuming 
R 0 = 10 kpc). This relatively high value for Q is in part due to the 
high mass density found by Innanen at R = R 0 and z=0 (0.149 I^/pc ) and 

in part to the approximations employed in deriving equation (20). All 
three curves exhibit the sharp rise in zjy2( R ) found by McGee and Milton 

beyond R 0 . These results are consistent with the following interpreta- 
tion: the increase in the half- thickness of the gas layer is due to 
a systematic decrease of p To at distances beyond R 0 , with Q remaining 
essentially constant. 

We have implicitly assumed here that the plane of maximum gas 
density is perfectly flat. 21-cm observations (Kerr 1957; Gum, Kerr, 
and Westerhout 1960) indicate, however, that the layer of interstellar 
hydrogen is warped. The points of maximum hydrogen density in the 
region R < R 0 define a plane to within 50 pc. However, for R > R 0 
the layer curves upward on one side of the galaxy and downward on the 
opposite side, the deviation from the principal plane reaching 1 kpc 
in the very outer rim of the galaxy. Theoretical discussions con- 
cerning the physical mechanism responsible for the distortion include 
interaction of the galactic halo with intergalactic gas (Kahn and 
Woltjer 1959), interaction of the galactic gas disk with the Magellanic 
clouds at or near their present distance (Avner and King 1967), and 
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interaction with the Magellanic Clouds during a single close passage 
(Habing and Visser 1967; Hunter and Toomre 1969). If the increase 
in half- thickness is causally related to the warpage, our derivation 
may not be applicable. 

The question naturally arises as to whether the increase in 
half- thickness of the gas layer is a general feature of extragalactic 
systems, as it should be if caused by a decrease in pto(R) at large R. 
Low resolution 21-cm studies of extragalactic system (Roberts 1968, 

1969) at present offer no answer to this question. We can only hope 
to use the interstellar dust of the external galaxy as a gas tracer. 
Further, the system must be observed almost exactly edge on so as to 
minimize projection effects which confuse the issue. One favorable 
such galaxy, NGC 5866, has been observed by Burbidge and Burbidge (1960). 
Its dust layer extends very far from the center of the galaxy and 
remains remarkably flat and coplaner within a certain distance of the 
galactic center. Beyond this distance, both the half- thickness and 
warpage increase considerably with increasing R. 

To pursue (ii), i.e. to construct a galactic mass model, we are 

faced with the problem of choosing appropriate values of Q(R). What 

makes this problem exceptionally difficult is that Q can change if 
2 2 

either <Vtz>> B o /Pgo> or p c-ro/Pgo vary with R. Because we quite 
obviously lack this information, Q is taken to be independent of R for 
purposes of simplicity, and is determined by the constraints that 
PTo = 0*15 I^/pc and Zjy 2 = 268 pc in the solar vicinity. Under these 
restrictions, Q becomes 10.25 km/sec. 
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The resultant galactic mass model is given in Table 2 along with 
Innanen' s galactic mass model for purposes of comparison. That p^o is 
independent of R over the range 4 kpc S R S 10 kpc is a result of the 
constancy of the observed layer half-width and the assumed constancy 
of Q. Our mass densities are systematically less than those found by 
Innanen over this region, the discrepancy becoming more substantial 
as R decreases. The agreement at R=R 0 is forced by our choice of Q. 
Uncertainties of at least a factor of 1.5 in Q introduce uncertainties 
greater than a factor of 2 in the mass model. To this must be added 
the uncertainties in the observed values of The importance 

of our mass model at R < Rq rests not in the exact numerical values 
derived for p To (R) (most certainly incorrect) but in demonstrating that 
present-day galactic mass models (e.g. Innanen 1966) cannot be 
reproduced by simply assuming that the gas and star disks obey hydro- 
static and Poisson considerations and that z^/2( R ^ an< ^ Q( R ) remain 
constant over the range 4 kpc ^ R s» 10 kpc. If the mass densities 
quoted by Innanen over this range are even near correct, it is clear 
that we must part with at least one of our assumptions. The restriction 
Q(R) = constant appears to be the weakest link in the argument. Indeed, 
if Q is determined by the conditions z^/2 = z l/2 observed and PTo = 

PTo innanen’ we find that Q = 24.8 km/ sec at R = 4 kpc, significantly 
greater than its value at R=R 0 . The suggested increase in Q(R) with 

decreasing R may be due to (a) larger gas turbulent velocities, possibly 
caused by an increase in the rate of star formation (Schmidt 1970) 
and or supernovae activity, (b) an increase in the term B 0 (R)/8rrpgo(R) > 
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or (c) an increase in the term P c _ ro (R)/p go (R) • It is interesting in 

this respect that 21-cm profile widths are observed to increase 

with decreasing R (Schmidt 1956), possibly indicating enhanced gas 

turbulence. The significant point to emerge from this discussion is 

simply that since zi/ 2(R) I s observed to remain essentially constant 

over the range 4 kpc S R S 10 kpc, then from equation (20) it follows 

that Q(R) °c p^R). If Pj 0 (R) increases with decreasing R, as seems 

reasonable, then Q must also increase. Finally, it should be emphasized 

that the above analysis should not be applied within 4 kpc of the 

galactic center because the nuclear component, which has been neglected, 

will begin to dominate the disk component in its gravitational effect 

2 2 

on the gas layer, and because 6 = Q /<v* z > may become i 1 so that the 
approximations employed in deriving equation (20) no longer are valid. 

IV. summary 

The distribution of gas density with distance above the galactic 
plane Pg(z) is computed in Section II on the assumption that the gas 
and star disks obey hydrostatic and Poisson considerations and that the 
gas turbulent, magnetic, and cosmic-ray pressures are in equipartition. 
Comparing the theoretical gas distribution with Schmidt's observations 
of Pg(z) at the galactic tangential points enables us to derive the 

rms z turbulent gas velocity dispersion. This dispersion 

is near 7 km/sec only if we accept the Oort mass limit and choose 
R 0 = 10 kpc . 

The same theoretical considerations are employed in Section III 
to derive an expression for the half-thickness of the galactic 
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gas disk. Employing the Innanen galactic mass model and the observations 

of McGee and Milton relating to the hydrogen layer half-width in the region 

4 kpc £ R ^ 14 kpc, it seems likely that (i) the increase in half-width 

for R > R 0 is due to a decrease in gas plus star density at the galactic 

plane with increasing R, and (ii) since the half-width remains essentially 

1/2 

constant when 4 kpc S R i 10 kpc, Q must scale as (Pg 0 + P* 0 ) and there- 
fore increase with decreasing R. 
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TABLE 1 


z TURBULENT GAS VELOCITY DISPERSION 


P s o + P *o 
(%/pc 3 ) 

R o 

(kpc) 

<vL >1/2 

(km/ sec) 

0.089 

8.2 

0 

0.089 

10.0 

3.00 

0.150 

8.2 

3.51 

0.150 

10.0 

6.82 

R 

TABLE 2 

GALACTIC MASS MODELS 
p-jo(Innanen) 

PT 0 (Ke liman) 

(kpc) 

(%/pc 3 ) 

(%/pc 3 ) 

4.0 

0.872 

0.149 

6.0 

0.470 

0.149 

8.2 

0.269 

0.149 

10.0 

0.149 

0.149 

12.0 

0.062 

0.0251 

14.0 

0.006 

0.010 

16.0 


0.0056 
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FIGURE CAPTIONS 

1. The distribution p (z)/p ( X for three values of Q, obtained 

O 

by solving the Poisson and two hydrostatic equilibrium equations. 

A smooth curve through Schmidt's observed data ( ) is shown. 

2. The distribution P e (z)/p eo ( ) obtained by solving the gas 

O o '" 7 

hydrostatic equilibrium equation with Oort's I^. A smooth curve 

through Schmidt's observed data ( ) is shown. 

3. The theoretically determined half-thickness of the gas layer 

( ) for three values of Q. A mean curve through the observations 

of McGee and Milton ( ) is shown. 
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